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TWO-DIMENSIONAL  CONVOLUTIONS,  CORRELATIONS,  AND 

FOURIER  TRANSFORMS  OF  COMBINATIONS  OF  WIGNER  DISTRIBUTION 
FUNCTIONS  AND  COMPLEX  AMBIGUITY  FUNCTIONS 

INTRODUCTION 

Over  the  years,  a  number  of  properties  of  integrals  of 
products  of  complex  ambiguity  functions  (CAFs)  or  products  of 
Wigner  distribution  functions  (WDFs)  have  been  derived,  such  as: 
the  volume  constraint  of  magnitude-squared  ambiguity  functions 
[1;  page  308],  the  positivity  of  the  convolution  of  any  two  WDFs 
[2;  (106)],  and  Moyal's  theorem  involving  the  volume  under  the 
square  of  a  WDF  [3].  Now,  it  appears  that  these  are  very  special 
cases  of  a  general  class  of  two-dimensional  Fourier  transforms  of 
combinations  of  CAFs  and  WDFs  with  delayed  or  time-reversed 
arguments . 

We  begin  by  deriving  a  general  one-dimensional  transform 
relation  involving  two  arbitrary  complex  waveforms  and  their 
Fourier  transforms.  An  application  of  this  relation  to  energy 
density  spectra  yields  three  alternative  expressions  for  the 
output  correlation  of  a  filtered  time  function.  This  general 
transform  relation  is  also  the  basic  tool  for  setting  up  the  two- 
dimensional  transforms  that  are  the  subject  of  succeeding 
sections.  The  extreme  generality  of  the  two-dimensional 
relations  allows  for  a  large  number  of  special  cases;  some  of 
these  are  pointed  out,  but  undoubtedly  there  are  additional  ones 
not  listed  here. 
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When  we  begin  our  two-dimensional  transform  investigation,  we 
do  not  immediately  specialize  to  WDFs  or  CAFs.  Rather,  we  first 
consider  a  set  of  four  general  functions,  each  of  two  variables, 
all  of  which  are  related  to  each  other  by  Fourier  transforms.  We 
show  that  two-dimensional  Fourier  transforms  of  products  of  pairs 
of  these  general  functions  are  all  equal  to  a  common  value, 
although  that  value  cannot  be  expressed  in  any  simple  closed 
form.  These  relations  are  derived  for  convolution  type 
operations  as  well  as  for  correlation  operations. 

When  we  make  a  specialization  of  these  results  to  waveforms, 
relatively  simple  closed  form  results,  in  terms  of  products  of 
WDFs  and  CAFs,  are  obtained  for  these  two-dimensional  transforms. 
And  when  the  arguments  of  these  relations  are  further  specialized 
in  value  (such  as  zero),  some  of  the  currently  known  relations 
involving  CAFs  and  WDFs  result. 

Extensions  of  these  results  to  time  contracted  or  expanded 
arguments  are  made  in  the  appendices.  Again,  specializations  to 
waveforms  yield  closed  form  results,  in  terms  of  products  of  WDFs 
and/or  CAFs. 
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ONE-DIMENSIONAL  TRANSFORM  RELATIONS 

Function  g(t)  is  an  arbitrary  complex  function  of  real 
argument  t,  which  will  be  thought  of  as  time.  Its  Fouriei 
transform  will  be  denoted  by  complex  function  G(f),  where 

G ( f )  -  Jdt  exp(-i2nft)  g(t)  .  (1) 

Integrals  without  limits  are  along  the  real  axis  and  over  the 
range  of  nonzero  integrand.  Argument  f  is  a  real  cyclic 
frequency,  not  a  radian  frequency.  The  inverse  Fourier  transform 
relation  to  (1)  is 

g(t)  -  Jdf  exp(fi2nft)  G(f)  .  (2) 

The  Fourier  transform  pair  in  (1)  and  (2)  will  be  denoted  by 

g(t)  *  G ( f )  .  (3) 

Similarly,  h(t)  and  H(f)  will  be  a  Fourier  transform  pair. 

TRANSFORM  OF  PRODUCT  OF  WAVEFORMS 

The  variables  v,a,@,/u,y  are  all  real  in  the  following.  A 
generalization  of  Parseval's  theorem  is  then  possible,  namely 

Jdt'  exp( -i 2 nvt '  )  g(at  +  |3t')  h*  ( //t+y  t '  )  -  exp  2 
x  Jdv'  exp^i2n\>'t(ay-@*y)j  G  (y  (V  +  )  h*(p(v'~  2^7) )  '  ^  > 
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where  it  is  presumed  that  3  *  0  and  y  #  0.  This  result  may  be 
derived  by  substituting  for  g  according  to  (2),  interchanging 
integrals,  and  using  (1)  for  Fourier  transform  pair  h(t)  «*  H(f). 
A  more  symmetric  form  for  relation  (4)  is  available,  if  desired: 

Jdt-  exp(-i2nvt'  )  g(e(t'+  j^))  h*(r(t'-  j^))  - 

-  Jdv  exp( +i2xv' t )  g(t(v.  jf-))  H*(e(v'-  .  (5) 

SPECIAL  CASES 

By  specializing  the  parameter  values  in  (4),  several 
interesting  and  useful  results  can  be  obtained.  For  example,  if 
we  take  y  -  8,  /u  -  -a,  then  we  obtain  a  combined  one-dimensional 
Fourier  transform  and  correlation: 

Jdt'  exp(-i2nvt')  g(0t'+at)  h*  (  f3t ' -at )  - 

-  Jd\> '  exp  (i2n\)'  t2a@j  g(sv'  +  H*^8v'-  .  (6) 

On  the  other  hand,  if  we  take  y  =  -8,  //  *  a  in  (4),  there 

follows  a  combined  one-dimensional  Fourier  transform  and 
convolution : 

Jdt'  exp(-i2nvt')  g(at+8t')  h*(at-8t')  - 

-  Jdv '  exp  ^i  2  nv '  t2  a8j  G  +8v'j  H*  -8v'j  .  (7) 
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Further  specialization  to  the  specific  numerical  values 
Y  ™  fi  m  1  *  ~/j  m  a  "  H  t  in  (6)  yields 

Jdt'  exp(-i2nvt')  g(t'+*st)  h*(t'-ijt)  = 

-  Jdv'  exp(+i2nv't)  G(v'+ijv)  H*(v'-^v)  .  (8) 

Alternatively,  the  choice  -y  «  8  =  H,  //  *  a  «  1  in  (7)  yields 

Jdt'  exp( -i2n\>t '  )  g(t+*?t')  h*(t-*st')  = 

-  Jdv'  exp( +i2  Jtv' t )  G(v+^v')  H*(v-^v')  .  (9) 

APPLICATION  TO  ENERGY  DENSITY  SPECTRA 

Case  1 .  Suppose  that  we  choose 

G( v )  -  |X(\>)|2  ,  H ( v )  -  | Y ( v ) | 2  ,  (10) 

which  are  the  energy  density  spectra  of  waveforms  x(t)  and  y(t), 
respectively.  Then  g(t)  -  ij»xx(t)  and  h(t)  -  ^yy(t),  where  *xx(t) 
is  the  auto-correlation  function  of  complex  waveform  x(t): 

*  (t)  '  fdu  x(t  +  u)  x*(u)  .  (11) 

XX  J 

The  use  of  (10)  and  (11)  in  (8)  yields 
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I ( t , v )  ■  Jdv'  exp(+i2nv't)  |X(v'+»{v)|2  |Y(v'-»iv)|2  - 

"  Jdt'  exp(-i2rivt' )  \|»xx(  t  '+»jt )  'I'yy  ( t ' -*st )  .  (12) 

The  last  term  in  (12)  is  identical  to  tyy(*st-t'). 

The  special  case  of  v  -  0  in  (12)  reduces  to 

Il{t,0)  "  Jdv'  exp(+i2nv't )  j X ( v ' ) | 2  |Y(v')|2  - 

“  Jdt'  *xx(t'+l5t)  ^yy  ( t f  —  1 )  •  (13) 

The  additional  restriction  to  t  -  0  becomes 

1^0,0)  -  Jdv'  |  X(  v'  )  |2  |  Y(  v'  )  j  2  *» 

-  jdt'  *XXU')  ^^y  (  t  '  )  .  (14) 


Case  2.  Here,  instead,  make  the  identifications 


G(  v )  -  X(  v )  Y(  v )  -  H (  v  )  . 


(15) 


?hen 


g(t)  -  C  (t)  ■  fdu  x(u)  y(t-u)  -  h(t)  , 

Xy  J 


(16) 


which  is  the  convolution  of  x(t)  and  y(t).  Substitution  of  (15) 
and  (16)  in  (8)  gives 
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I2(t,v)  ■  Jdv'  exp( +i2nv' t )  X(v'+%v)  Y(v'+4v)  X*(v'-»{\>)  Y*(v'-Sjv) 

«  fdt'  exp( -i2nvt '  )  C  (t'+Ht)  C*  (t'-Jjt)  .  (17) 

J  xy  xy 

Setting  v  to  zero  yields 

l2(t,0)  -  Jdv'  exp(+i2rtv't)  | X ( v ' ) | 2  |Y(v')|2  « 

-  Jdt'  C  (t'+%t)  C*y  ( t '  — *5 1 )  .  (18) 

Finally,  also  setting  t  equal  to  zero, 

t2(0,0)  -  Jdv'  |X(v')|2  | Y ( v ' ) | 2  -  Jdt'  |Cxy(t')|2  .  (19) 

Case  3.  Now  identify 

G( v )  -  X(v)  Y* ( v )  -  H ( v )  .  (20) 

Then 

g(t)  =  ( t )  *  Jdu  x(u  +  t)  y* ( u )  -  h ( t )  ,  (21) 

which  is  the  cross-correlation  of  x(t)  and  y(t).  The  use  of  (20) 
and  (21)  in  (8)  leads  to 

Ij(t,v)  ■  Jdv'  exp(  +  i2n\>' t )  X(v'+Jjv)  Y*(v'+*{v)  X*(v'-Hv)  Y(v'-^v) 

-  Jdt'  exp( -i2rcvt ' )  ^xy(t'+ht)  **  ( t ' -%t )  .  (22) 

The  result  of  setting  v  to  zero  is 
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I3( t,0)  -  Jdv'  exp(+i2nv't)  |X(v')|2  |Y(v')|2  ■ 

-  Jdt '  *xy(t'+««t)  **y(t'-Ht)  .  (23) 

When  t  is  also  set  equal  to  zero,  (23)  reduces  to 

1 3  (  0 , 0 )  -  Jdv'  |X( v' ) | 2  |  Y(  v'  )  |  2  -  Jdt'  M'xytt')!2  .  (24) 

It  should  be  observed  that  the  upper  lines  of  (13),  (18),  and 
(23)  are  identical  to  each  other;  that  is, 

l1(t,C)  -  l2(t,0)  -  l3(t,0)  .  (25) 

Therefore,  the  lower  lines  of  (13),  (18),  and  (23)  furnish  three 
equal  alternative  expressions  involving  autocorrelations, 
convolutions,  or  cross-correlations,  respectively. 

There  are  many  other  possibilities  for  identifications  of  G 
and  H  in  (8),  besides  (10),  (15),  and  (20).  For  example,  we 
could  take 

G( v )  -  | X( v) | 2  Y(v)  ,  H( v )  -  Y ( v )  .  (26) 

However,  it  may  be  shown  that  this  choice  leads  identically  to 
result  (13)  when  v  is  set  to  zero;  so  not  all  selections  yield 
new  relations.  Additional  convolution  type  relations  may  be 
obtained  if  (9)  is  used  instead  of  (8). 
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GENERAL  TWO-DIMENSIONAL  TRANSFORM  RELATIONS 

In  this  section,  we  will  consider  a  set  of  four  general 
functions,  each  of  two  variables,  which  are  related  to  each  other 
by  Fourier  transforms.  These  four  functions  are  indicated  in 
figure  1,  where  a  two-headed  arrow  denotes  a  Fourier  transform 
relationship.  These  functions  are,  for  the  moment,  arbitrary 
complex  functions  of  two  variables;  they  are  not  necessarily 
Wigner  distribution  functions  or  complex  ambiguity  functions. 

R(  t ,  x )  <- - ->  X(v,t) 

i  i 

W(  t ,  f  )  < - *(  v,  f  ) 

Figure  1.  General  Two-Dimensional  Functions 

The  paired  transform  variables,  here  and  for  the  rest  of  the 
report,  are  t  «*  v  and  x  **  f .  The  detailed  Fourier  transform 
interrelationships  between  the  four  functions  in  figure  1  are 


X(v,t)  -  Jdt  exp(-i2nvt)  R(t,x)  ,  (27) 

R ( t , x )  »  Jdv  exp( + i2nvt )  x(v,x)  ,  (28) 

W(t,f)  =  Jdx  exp(-i2nfr)  R(t,x)  ,  (29) 

R(t,x)  =  Jdf  exp(+i2nfx)  W(t,f)  ,  (30) 


9 


TR  8759 


♦<v,f)  - 

Jdt 

exp( -i2nvt ) 

W(t,f)  , 

(31) 

W(t,f)  - 

Jdv 

exp( +i2revt ) 

*(v,f)  , 

(32) 

*(v,f)  - 

JdT 

exp( -i2nf x) 

X(v,x)  , 

(33) 

X(v,x)  - 

Jd£ 

exp( +i2nfx) 

♦(v,f)  . 

(34) 

A  double  Fourier  transform  relationship  exists  between  K  and  #, 
as  well  as  between  W  and  x* 

TWO-DIMENSIONAL  CONVOLUTIONS 

We  repeat  (9)  here,  but  with  a  change  of  variables  t  -*  t  and 
v  -*  f : 

Jdx'  exp( -i2nf x' )  g(x+Hxr)  h*(x-*jx')  - 

-  Jdf '  exp(+i2nf ' x)  G(f+Hf')  H*(f-«sf')  .  (35) 

Let  X^  and  X2  be  two  different  functions  of  the  type  indicated  in 
figure  1,  and  consider  (35)  with  the  assignments 

g(x)  -  X1(va,x)  ,  h(x)  -  X2<vb,x)  .  (36) 

The  corresponding  Fourier  transform  pairs  for  (36)  are 

G(f)  -  VV*)  '  H(f)  "  *2(vb'f)  '  (37) 

upon  use  of  (33).  There  follows,  from  (35), 
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Jdx*  exp( -i2nf x'  )  va'T+l«T' >  Xj  ( ^b*  )  - 

■  Jdf'  exp(  +i2itf' t)  (  va ,  f+*$f  •  )  **(vb,f-4f')  .  (38) 

See  appendix  A  for  the  most  general  result  of  this  form. 

If  we  now  let  \>a  -  v+^v'  and  \>b  -  v-*sV  in  (38),  then  an 
additional  Fourier  transform  on  v'  yields  the  middle  two  lines 
in  (39)  below.  More  generally,  in  a  similar  fashion  to  that  used 
above,  we  find  that  the  combined  two-dimensional  convolution  and 
Fourier  transform  can  be  expressed  in  four  equivalent  forms: 

I(\>,f,t,x)  ■  (39) 


-  JJdt'dT' 


dv'  dx ' 


-JJ 

-Jf 


dv'df ' 


dt'df ' 


exp(-i2nvt'-i2nfx'  )  R^  ( t+*st '  ,  x+>jx'  )  ( t-*jt '  ,  x-*sx'  ) 
exp( +i2nv' t-i2nf x' )  Xj^  (  v+*iv'  ,  x+^x'  )  X^  (  v-*sv'  ,  x-*sx'  ) 
exp(  +i2nv'  t+i2nf  ’  x)  ( v+*sv' ,  f+^f '  )  *2  (  v-*}v' ,  f-»sf '  ) 
exp(-i2nvt'4i2nf rx)  ( t+Ht '  ,  f ( t— *s t ' ,  f — *5 f '  ) 


Alternative  forms  of  (39)  are  available;  for  example,  the 
last  line  can  be  written  in  the  more  typical  convolution  form 

JJdt'df'  exp( -i2nvt ' +i 2nf ' x )  W^(t',f')  ( t-t ' , f-f ' )  - 

■  >*  exp(-invt+infx)  I  (*{v,iif  ,*st,*sx)  .  (40) 
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TWO-DIMENSIONAL  CORRELATIONS 


Here,  we  use  (8)  with  identifications 


g(t)  -  R1(t,xa) 
G(v)  -  Xx(v,xa) 

Then  there  follows  immediately 


h(t)  -  R2(t,xb)  , 
H(  v )  *=  X2(v,xb)  • 


exp(-i2revt')  R1{t'+ht,T  ) 
x  a 

exp<+i2re\>'t)  X1( v'+*iv,xa) 


R2^  t,-**t#  Tfa)  - 
X2(  v'->jv,Tb)  . 


(41) 


(42) 


Now  let  xfl  -  t'+*st  and  xb  -  x'-hx,  and  Fourier  transform  on 
x' .  The  result  is  the  first  two  relations,  given  below,  of  four 
equivclent  forms  of  the  combined  two-dimensional  correlation  and 
Fourier  transform 


■  Jjdt'dT 

-  Jjd',’dT 

- 

-  JJdt’df 


r 


f 


F 


F 


J(v,f ,t,T) 


exp(  -i2xvt' -i2nf  x' )  R^  ( t'+»jt  ,x'+^x) 


exp(+i2nv'  t— i 2 nf  "c ' )  Xj^  ( vr  +4v,  x'  4-Jjt) 


exp( +i2nv' t+i2nf r  x)  ^  ( v' +»jv,  f ' +%f ) 


exp(  -i2?ivt *  +  i2 nf '  x)  ( t ' 4-Jjt ,  f  r+»sf ) 


(43) 

R*(  t'-ijt,x'-»sx)  - 
X2<  x'-*jx)  - 

«2(v'-Jjv,f'-»jf )  - 
W*(t'-»5t,f'->jf )  . 


Alternative  forms 
line  can  be  expressed 


to  (43)  are 
in  the  more 


possible;  for  example,  the 
typical  correlation  form 


last 
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JJdt'df'  exp( -i2nvt '+i2nf ' t)  W1(t',£')  ( t ' -t , f ' -f )  - 

-  exp(-invt+infT)  J(v,f,t,t)  .  (44) 

MIXED  RELATIONS 

The  results  in  (39)  and  (43)  all  involve  two  W(t,f) 
functions,  or  two  X(v,t)  functions,  etc.  However,  it  is  possible 
to  obtain  relations  which  involve,  for  example,  one  W(t,f) 
function  and  one  x(v,t)  function.  As  an  illustrative  example, 
consider  (9)  with  g(t)  -  W1(t,fa)  and  h(t)  «  x2  ^  )  •  Then, 

from  figure  1,  G(v)  -  anc*  “  *2^fb'v*'  giving 

Jdt '  exp(  -i2rtvt '  )  ( t+«jt '  ,  f  & )  X^  (  t-«jt '  )  - 

-  Jdv'  exp(+i2nv't)  *^ ( v+Hv' , f g f^,  v-Jjv' )  .  (45) 

If  we  now  let  f  ■  f+*jf'  and  f.  -  f-*s£',  and  perform  a 

3  D 

Fourier  transform  on  f',  there  follows  immediately 

JJdt'df'  exp( -i2nvt '  +  i2rcf  'x)  ( t+*jt '  ,  f -»-*4  f '  )  x^  (  f — *5  f  ’  ,  t-*jt '  )  - 

(46) 

«  JJdv'df'  exp(  +i2rev'  t+i2ref '  x)  ( v+*jv' ,  f +%f '  )  ♦  ]£(  f-**f '  ,  v-*jv'  )  . 

Thus,  a  combined  two-dimensional  convolution  and  Fourier 
transform  of  a  W(t,f)  function  and  a  X(v,t)  function  can  be 
expressed  in  terms  of  two  *(v,f)  functions.  (Strictly,  some  of 
the  arguments  are  reversed,  as  seen  in  (46).) 
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If,  instead,  we  use  (8)  with  g(t)  and  h(t)  assigned  as  above, 
then  we  obtain 

Jdt'  exp( -i2nvt ' )  W1(t,+%t,fa)  *2*  fb' - 

-  Jdv'  exp(+i2nv't)  *1  ( vr+«jv,  fft )  ♦*  (  v'-»jv)  .  (47) 

Letting  «  f'+*jf,  «  f'-Hf,  and  performing  an  additional 

Fourier  transform  on  f',  there  follows 

[Jdt'df'  exp(-i2nvt'  +  i2xf 't)  t'+«jt ,  f  '+>jf )  xj<f  r-*|f  #t'-»it)  - 

(48) 

-  JJdv'df'  exp(  +i2nv'  t+i2nf  •  x)  ♦1(  v'+%v,  f  r+*jf )  f  '-Hf ,  v'-*jv)  . 

Here,  a  combined  two-dimensional  correlation  and  Fourier 
transform  of  a  W(t,f)  function  and  a  X(v,t)  function  can  be 
expressed  in  terms  of  two  *(v,f)  functions.  (Again,  some 
arguments  are  reversed  or  replaced.  However,  the  first  argument 
in  a  X  function  is  always  a  frequency  variable,  while  the  second 
argument  is  always  a  time  variable;  similar  restrictions  hold  for 
the  remaining  functions  R,  W,  ♦  in  figure  1.) 
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SPECIALIZATION  TO  WAVEFORMS 


In  the  previous  section,  the  functions  R,  W,  X,  *  were 
arbitrary,  except  that  they  were  related  by  Fourier  transforms 
according  to  figure  1.  Here,  we  will  specialize  their  forms, 
thereby  enabling  more  explicit  relations  for  their  two- 
dimensional  convolutions  and  correlations. 

For  arbitrary  complex  waveforms  a(t),  b(t),  c(t),  d(t),  let 

R^t,!:)  =  a(t+Hx)  b*(t-*5T)  *  Rab(t,t)  ,  (49) 

R2(t,T)  *  C(t+^T)  d*  (  t-^T)  a  Rcd(t,T)  .  (50) 


These  are  known  as 
Thus,  Rab(t,x)  is 
waveforms  a  and  b, 
(or  delay)  time  x. 
follows 


(cross)  temporal  correlation  functions  (TCFs). 
the  "instantaneous”  cross-correlation  between 
corresponding  to  center  time  t  and  separation 
Then,  from  (31)  and  (29),  or  [4;  (35)],  there 


*1(v,f)  «  *ab( v, f )  -  JJdt  dx  exp( -i2nvt-i2nf x)  Rab(t,r) 


=  A(  f  +  ^v)  B  (  f->{V )  , 


(51) 


*2(v,f)  -  *c<J(v,f)  -  C(f+*sv)  D*  (  f-*«v )  .  (52) 

These  functions  are  known  as  (cross)  spectral  correlation 
functions  (SCFs).  (in  (4),  the  notation  A(v,f)  was  used  for  this 
function;  however,  A(f)  will  be  used  here  for  the  Fourier 
transform  of  waveform  a(t).)  The  SCF  corresponds  to  center 
frequency  f  and  separation  (or  shift)  frequency  v. 
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The  Fourier  transform  relationships  in  figure  1  and  equations 
(27)  -  (34)  still  hold  true,  but  now  are  specialized  to  the 
waveform  cases  above.  Specifically,  figure  2  illustrates  the 
four  two-dimensional  functions  for  waveforms  a(t)  and  b(t),  where 
now  **  Waj3  is  a  cross  Wigner  distribution  function  (WDF)  and 
Xf  ■  Xab  is  a  cross  complex  ambiguity  function  (CAF). 


TCF  Rab(t,x)  <- - ->  Xab(\>,T)  CAF 


WDF  Wab(t'f)  * - *ab(v,f)  SCF 


Figure  2.  Two-Dimensional  Functions  for  Waveforms 


The  detailed  Fourier  transform  interrelationships  are  now 


Xab(v,T*  "  Jdt  exP<-i2rtvt)  Rab(fc»T)  ,  (53) 

Rab*t,T)  "  Jdv  exP(  +  i2l*vt)  Xab(v,x)  ,  (54) 

Wab(t,f)  *  Jdr  exp(-i2nft)  Rab(t»T)  ,  (55) 

Rab*t,T*  "  Jdf  exP^  +i2nf t)  Wab(t,f)  ,  (56) 

*ab(v,f)  "  Jdt  exP(-i2nvt)  Wafa(t,f)  ,  (57) 

WabU,f)  -  Jdv  exp(  +i  2  nvt )  *ab(v'f)  *  (58) 
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Jdx  exp(-i2nfx) 

Xab(V.T)  , 

(59) 

Xab<v,t)  - 

Jdf  exp(+i2itfx) 

4ab<'''£>  • 

(60) 

The  function  W__(t,f),  for  example,  is  an  auto  WDF ,  since  it 
involves  only  one( waveform,  a(t).  We  will  frequently  drop  the 
terminology  auto  and  cross,  when  possible  without  confusion,  and 
let  the  notation  indicate  the  particular  case. 

It  will  .be  found  advantageous  for  future  purposes  to  define  a 

t 

scaled  and  contracted  WDF  according  to 

Wab(t,f)  -  *«Wab(  »st ,  »sf  )  .  (61) 

GENERAL  CROSS  PROPERTIES 

Due  to  the  restriction  of  form  taken  on  by  the  TCF  in  (49) 
and  the  SCF  in  (51),  the  four  functions  in  figure  2  obey  some 
symmetry  rules;  they  are 


Rab<t'-T> 

‘  Bba(t'T)  ' 

• 

♦ab<-v'£) 

■  ♦ba(v’£)  ■ 

• 

Xab<-'’--'> 

’  Xba(  v»  x)  , 

•  "ba(t'£)  • 

(62) 
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AUTO  PROPERTIES 

When  waveform  b(t)  *  a(t),  some  specializations  follow: 


aa(t'-T) 

*  Raal,:'T)  - 

-  Cl','£)  • 

:aa(-','-T) 

*  Xaa(v-T)  ' 

aa't'f'  - 

real  for  all  t,  f,  a(t), 

(83) 

with  the  only  significant  specialization  being  the  realness  of 
WDF  W  (t,f).  Waveform  a(t)  can  still  be  complex. 

O  u 

SOME  SPECIAL  CASES 

The  ordinary  cross-correlation  of  two  waveforms  a(t)  and  b(t) 
is  a  special  case  of  a  CAF: 

*ab(T)  ■  Jdt  a(t)  b*(t-x)  -  xab*°'T*  *  ^64) 

The  ordinary  cross-spectrum  is  then  a  special  case  of  an  SCF: 

Yab(f)  ■  Jdr  exp(  -i2nf  t)  \j/ab(x)  «  ^ab*0'1^  "  B(f)  •  (55) 

The  autospectrum  is  then  simply 

Y(f )  -  *  (0,f  )  -  j  A  (  f )  |  2  ,  (66) 

<3  cl  da 

which  is  always  nonnegative. 

The  ordinary  convolution  of  two  waveforms  a(t)  and  b(t)  is  a 
special  case  of  a  WDF: 

JdT  a(T)  b*  ( t-r)  -  *iWab(«*t,0)  -  Wab(t,0)  .  (67) 
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REAL  WAVEFORM  a(t) 

In  addition,  if  waveform  a(t)  is  real,  the  following  (auto) 
properties  hold  true: 


The  situation 
below. 


aa^'-O 

-  R__(t,T)  and  R=„  is  real  , 

ad  aa 

aa^'-f) 

-  «aa<"'£>  ’ 

aa(v'-T) 

-  xaa<','r>  ' 

-  “aa*1'11  ' 

(68) 

for  a  real  waveform  a(t)  is  summarized  in  figure  3 


Figure  3.  Symmetry  Properties  for  Real  Waveform  a(t) 
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MIRROR-IMAGE  RELATIONS 

For  general  complex  waveforms  a(t)  and  b(t),  define 
mirror-image  functions 

a(  t )  -  a(-t )  ,  b(t)  -  b(-t)  .  (69) 

Then  it  follows  directly  that  the  voltage  density  spectrum  of 
mirror-image  a(t)  is 

A(f)  ■  Jdt  exp(-i2nft)  a(t)  -  A(-f)  ,  (70) 

which  is  the  mirror-image  of  A(f).  Also,  there  follows 


"  Rab*t,1c*  ' 

-  *ab<v'f)  ' 

Xab(-v,-T) 

-  X?t,<v,t)  , 

“  wab(t'f)  ' 

Thus,  the  mirror-image  property  for  A(f)  carries  over  into  all 
the  two-dimensional  domains,  such  as  the  WDF  and  CAF,  as  well. 
There  is  no  significant  simplification  for  b(t)  -  a(t),  except 
for  the  realness  of  W  (t,f),  as  before. 

99 

Use  of  mirror-image  definition  (69)  allows  for  an  interesting 
connection  between  WDFs  and  CAFs.  First,  substituting  (49)  into 
(53)  and  (55),  we  have  cross  CAF 
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xab*v'T*  "  Jdt  exP(~i2nvt)  a(t+»jx)  b*(t-*5X)  - 

-  Jdf  exp( +i2nf x)  A(f+»*v)  B*(f-Jjv)  ■  XA0(v»x)  (72) 

and  cross  WDF 

Wab(t»f)  ”  Jdx  exp(-i2nfx)  a(t+*sx)  b*(t-*$x)  « 

■  Jdv  exp(+i2nvt)  A(f+%v)  B*(f-»{v)  ■  wAB(t»f)  •  (73) 

Reference  to  (69)  now  immediately  reveals  that 

wab<t,f)  -  2Xab(2£,2t)  (74) 

or 

xab(v,T)  “  JSwab(,5T',{V)  "  Wftb(  x,  v )  .  (75) 

Here,  we  also  used  (61).  That  is,  the  WDF  of  two  waveforms 
a  and  b  is  proportional  to  the  CAF  of  waveforms  a  and  b,  the 
mirror-image  of  b. 

Finally,  since 

B* ( f )  *  b* ( -t )  -  b*  ( t )  ,  (76) 

then,  using  ( 72 ) , 

X  «  [df  exp(i2rtfx)  A(f+*jv)  B(f~4v)  * 

AB  J 

-  X  *(>>,x)  *  W  *(lJX,»sv)  -  W  *(x,v)  .  (77) 

ab  ab  ~ab 
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TWO-DIMENSIONAL  TRANSFORM  RELATIONS  FOR  WAVEFORMS 

In  an  earlier  section,  general  two-dimensional  transform 
relations  were  derived  between  sets  of  four  functions  related  by 
Fourier  transforms;  see  figure  1  and  (39)  and  (43).  Here,  we 
will  utilize  the  particular  forms  considered  in  the  previous 
section  for  waveforms  (see  figure  2)  and  will  derive  closed  forms 
for  I  and  J  in  (39)  and  (43),  respectively. 


TWO-DIMENSIONAL  CONVOLUTIONS 

If  we  substitute  (49)  and  (50)  in  the  top  relation  in  (39), 
there  follows 


I(v,f,t,x)  -  JJdt'dx'  exp(-i2nvt'-i2nfx'  )  a(  t+^t  '+*sx+>*x'  )  x 
x  b*(  t+*st ' -%x-(<x'  )  x  c*  ( t-*jt  > )  d(  t-*jt  '-J?x+%x'  )  .  (78) 


Now  let 


u  -  *{t'+*<x',  v  -  ijt'-Hrr';  u+v  «  t',  2(u-v)  *  x'  . 


(79) 


Since  the  Jacobian  of  this  transformation  is  4,  (78)  becomes 

!(>,,£, t,r)  -  ijjdu  dv  exp  ^-i2it\>(  u+v  )-i2nf  2  { u-v )  j  x 
x  a(t+Hx+u)  b*(t->jx+v)  c*(t+4x-u)  d(t-*jx-v)  - 


-  Jdu'  exp(-i2nu'  (  f+*jv) )  a(  t+*jx+»ju' )  c*  ( t+^x-^u' )  x 
x  Jdv'  exp( +i2nv'( f-4v) )  b* ( t-Hx+Sv ' )  d( t-»jx-**v ' )  - 
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^ac(  t+Sx,  f+%v)  wbd(  f-**v)  . 


(80) 


That  is,  all  the  following  quantities  are  equal 

I(v,f,t,r)  - 


■  JT 

■If 

-u 

-n 


dt'dt' 

exp( -i2nvt '-i2nf t' ) 

F‘ab  ( t+*5t '  ,  t+Jjt '  ) 

R*d( t-«ftr  ,x-Hx'  )- 

i 

dv'dr' 

exp( +i2nv' t-i2nf r' ) 

Xab(v+»jv'  ,x+Hr’  ) 

X*d<V-*V  )- 

dv'dt ' 

exp( +i2  ftv' t+i2nf '  x ) 

*ab{v+1*v'  'f+1*f ' ] 

^d(v~*lv''f~‘s£')“ 

dt'df » 

exp( -i2nvt '  +  i2nf '  x) 

wab(  t+*st ' » f+%f ' ) 

w*d(  t-*jt '  ,  f-*jf  f  )- 

-  W9C(  t+‘sT,f4‘*v) 

(81) 

All  four  double-integrals  in  (81)  can  be  expressed  as  a  product 
of  the  same  two  one-dimensional  integrals,  which  are  cross  WDFs. 
This  reduction  is  only  possible  when  the  two-dimensicnal 
functions,  like  Wab  and  Xab,  are  WDFs  and  CAFs,  respectively. 

The  transformations  in  (81)  are  combined  two-dimensional  Fourier 
transforms  and  convolutions  of  TCFs,  CAFs,  SCFs,  or  WDFs. 

By  use  of  (74),  an  alternative  expression  for  the  end  result 
in  (81)  is 

I(v,f,t,r)  -  4  Xac ( 2f +v , 2 t+T )  X^d( 2f-v,2t-x)  ,  (82) 

in  terms  of  mirror-image  functions;  see  (69).  Also,  a  more 
typical  convolution  form  for  (31),  for  example,  is  (using  (61)) 
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JJdu  dv  exp( -i2n\>u+i2  nvx)  ( u ,  v )  W*d(  t-u,  f-v )  « 


-  exp(-invt+infx)  Wac(  t+>jx,  f+»jv)  wfad(  t-»sx,  f-»iv)  . 


bd' 


(83) 


TWO-DIMENSIONAL  CORRELATIONS 


In  an  identical  fashion  to 

becomes 

that  used  above, 

result  (43) 

J(  v, 

f,t,x)  - 

•  I. 

Jdt'dx'  exp( -i2n\»t ' -i2nf x’ ) 

Rab{  t'+l*t'T,+,5x) 

R*d(  t'-*jt,x'-»ix)- 

*  I 

Jdv'dx'  exp( +i2nv' t-i2nf x' ) 

Xab(  V'+,*V'T'+,5T) 

X*d(  v'-»*v,x'-»ix)- 

■I 

Jdv'df'  exp( +i2nv' t+12rtf ' x) 

*ab(  v'+ljV'f  '+**f  > 

♦  *d(  f 

■  I 

jdt'df'  exp( -i2nvt ' +i2 nf ' x) 

wab(t'+,jt'f'+**£) 

-  Xgc(  f  +  *JV,t  +  »JT) 

Xbd(f-‘iV,t-**x)  • 

(84) 

All 

these  double  integrals  in 

(84)  are  equal  to 

a  product  of  two 

cross  CAFs.  Again,  this  only  holds  for  the  special  forms  of  the 
two-dimensional  functions,  like  Wab  and  which  are  WDFs  and 

CAFs,  respectively.  The  transformations  in  (84)  are  combined 
two-dimensional  Fourier  transforms  and  correlations  of  TCFs, 
CAFS,  SCFS,  Or  WDFS. 

By  use  of  (75),  an  alternative  expression  for  the  end  result 
in  (84)  is 
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J(v,ff,t,T)  -  WftC(t+>jT,f-»-»*v)  Whd(t-»iT,f-»jv)  , 


bd 


(85) 


in  terms  of  mirror-image  functions.  Also,  a  more  typical 
correlation  form  for  (84)  is,  for  example, 

JJdu  dv  exp( -i2nvu+i2nvT)  Wftb(u,v)  W*d(u-t,v-f )  ■ 


-  exp( -i nvt+i nf  x)  Xftc<  f+hv,  t+Jjx)  Xbd(  f-^v,  t-»jT)  . 


(86) 


A  MIXED  RELATION 


As  an  example 
M1(t,f ) 


in  this  category,  if  we  take  (46)  with 
-  Mab(t,f)  ,  X2(v,t)  -  Xcd(2v,2t)  , 


then 


(87) 


*2_(\>,f)  *  #ab(v,f)  "  A(f+>sv)  B*(f-«jv)  , 

♦2(v,f)  -  H  *cd(2v,»«f)  -  >s  C(»sf-fv)  D*  ( *jf-v)  .  (88) 

Substitution  of  these  results  in  (46)  yields 

JJdt'df'  exp(-i2nvt'+i2xf 'x)  wab(  t+»*t ' , f ' )  X*d(2f-f '  ,2t-t' )  - 

-  Wflc(  t+HT,f+*iv)  Xbd(2f-v,2t-x)  .  (89) 

This  mixed  relation  is  a  two-dimensional  Fourier  transform  and 
convolution,  involving  a  WDF  and  a  CAF,  expressible  in  closed 
form  as  a  product  of  another  WDF  and  CAF. 
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SPECIAL  CASES 


The  two-dimensional  transform  results  in  (81)  and  (84)  in  the 
previous  section  involve  four  argunw ,.t ,  namely  v,f,t,x,  and  four 
functions,  a(t),  b(t),  c(t),  d(t).  Thti.  me  generality 

allows  for  numerous  special  cases  upon  seloc  •  >  >.  r*  c  f  the 
arguments  and/or  the  functions.  We  consider  cc-ne  .if  these 
possibilities,  but  are  aware  that  this  list  coulo  be  considerably 
augmented. 


Case  1.  As  an  example  of  the  generality  of  these  results, 
consider  in  (84)  the  particular  selection 


v-f-t-r-0,  c(t)  »  a(t),  d(t)  -  b(t). 


(90) 


There  follows  immediately  the  "volume  constraint" 


JJdv'dx'  Xab(v',T')  2  -  JJdt'df'  Wab(t',f 


')  2 


-  Xaa(0,0)  Xbb(0,0)  -  Jdt  | a ( t )  |  2  Jdt  | b ( t ) | 2  .  (91) 

Case  2.  in  (84),  take  v  »  x  «  0,  b(t)  -  a(t),  d(t)  -  c(t). 
Then  there  follows,  upon  use  of  (85), 


Jf 


dv'dx'  exp( +i2nv' t-i2nf x' )  x__(v',x')  x„„(vr,x')  - 

a  &  C  C 


dt'df '  w 


aa 


-  X 


ac 


(t'+Ht,f'+ijf )  w  (t'-»it,f»-«jf)  - 

C  w 

(f ,t) I2  -  |wac(t,f )  2  , 


(92) 
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which  is  nonnegative  real  for  all  f,  t,  a(t),  c(t).  Thus,  the 
two-dimensional  correlation  of  two  auto  WDFs  is  nonnegative. 

An  alternative  form  of  (92)  is 


fjdu  dv  Waa(u,v)  wcc(u-t,v~f)  - 


(93) 


Further  specialization  to  t  ■  f  ■  0  yields 

JJdu  dv  Waa(u,v)  Wcc(u,v)  -  j  Xac( 0,0) | 2  -  |jdt  a(t)  c*(t)|2,  (94) 
which  yields  Moyal's  result  (3)  for  c(t)  -  a(t),  namely 


JJdt  df  W2a(t,f)  -  [Jdt  | a ( t ) | 2] 2  . 


(95) 


Case  3.  In  (81),  take  v  -  x  -  0,  b(t)  -  a(t),  d(t)  -  c(t). 
We  then  get  the  "smoothing  result" 


JJdt'df'  Waa(t+»jt',f+»jf' )  Wcc(  t-»st'  ,f-»jf '  )  - 


-  wac(t'f> 


jjdx'  exp(-i2nfx' )  a(t+»jx')  c*(t-Hx')  2  >  0  (96) 


for  all  t,  f,  a(t),  c(t).  An  alternative  form  is 


JJdu  dv  Waa(u,v)  Wcc(t-u,f-v)  - 


*5Wac(‘st'J5f)  2  "  Wac(t,f )  2 


jjdx'  exp( -i2nf x' )  a(x')  c*(t-x')|2  . 


(97) 


That  is,  the  two-dimensional  convolution  of  two  auto  WDFs  is 
never  negative  (just  as  for  the  correlation  in  (92)). 
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Case  4 .  Using  (62),  the  same  basic  end  result  is  obtained 
from  (81)  for  the  following  double  integral  involving  CAFs: 


JJdv'dr'  exp(+i2nv' t-i2nfT' )  Xaa(*iv' ,%t' )  Xcc(iiv'  ,»jt'  ) 


wac(t'f) 


(98) 


This  right-hand  side  is  nonnegative  real  for  all  t,  f,  a(t), 
c(t).  An  alternative  form  is,  upon  use  of  (61), 


fj 


dv  dT  exp( +i2nvt-i2nf r) 


xaa(v'T>  Xcc(v'T) 


(99) 


Case  5.  Consider  (81)  with  c(t)  -  a(t),  d(t)  -  b(t).  Then 
the  right-hand  side  of  (81)  is  always  real.  For  example,  we  have 


JJdv'dt'  exp(  +i2nv'  t-i 2nf x' )  xat)(  v+%v' ,  t+*st' )  x*b(  v-ijv'  ,  )■■ 


«*  JJdt'df'  exp(-i2nvt'+i2nf '  t)  Wab(  t+Ht '  ,  f +*jf '  )  W*b(  t-*jt'  ,  f — f '  )■ 


"  Waa^ t+ljT'  Wbb(t-»sT,f-Hv)  .  (100) 

This  is  real  for  all  t,  t,  f,  v,  a(t),  b(t),  although  it  could  go 
negative . 

Case  6.  From  (81),  with  v  «  t  -  0,  there  follows 

JJdt'df'  Wab(t+ht' ,f+»if '  )  w*d(t— »st'  ,f-»jf '  )  - 

-  wac(t,f)  wbd(t,f)  ,  (101) 
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or,  with  the  help  of  (61)  and  (75),  alternative  form 

JJdu  dv  Wafe(u,v)  W*d(t-u,f-v)  - 

-  ?ac(t,f)  W*d(t,f)  -  Xac(f,t)  Xbd(f,t)  .  (102) 

Furthermore,  if  we  set  c(t)  -  a(t),  d(t)  -  b(t),  we  obtain 

JJdu  dv  Wab(u,v)  W* b ( t— u , f — v )  - 

-  ^bb(t'£)  •  V'£,t>  Xbb(£'U  •  (103> 

Thus,  the  two-dimensional  convolution  of  a  complex  cross  WDF  with 
itself  is  always  real,  but  could  go  negative. 

Case  7.  From  (83)  and  (84),  with  v  -  x  »  0,  there  follows 

JJdt'df'  wab(t'+«jt,f  '+Jjf )  W*d  ( t '  — t ,  f r  -  Jj  t )  - 

«  JJdu  dv  Wab(u,v)  W*d ( u-t , v-f )  - 
-  JJdv'dt'  exp( +i2nv' t-i2sf x' )  xab(v',r')  X*d(v',T')  - 

"  xac(£'U  xbd(f'fc)  ■  *ac(t'f)  *bd(t'£)  *  (104) 

The  two-dimensional  correlation  of  two  cross  WDFs  is  a  product  of 
two  cross  CAFs. 
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Case  8.  If  we  now  set  c(t)  -  a(t)  and  d(t)  -  b(t)  in  (104), 
we  obtain 

JJdu  dv  wab<u*v)  W*b(u-t,v-f)  “ 

*  JJt*v'<^‘r'  exp(  +i2  nv'  t-i2nf  x'  )  Xab(v',x')  2  - 

“  X*b(f,t)  -  Wa§(t,f)  w;b(t,f)  .  (105) 

The  two-dimensional  correlation  of  a  cross  WDF  with  itself  is  a 
product  of  two  auto  CAFs. 

Case  9.  From  (84),  with  t  -  f  -  0,  c(t)  -  a(t),  d(t)  -  b(t), 
and  with  the  help  of  (63),  we  find 

JJdt'df'  exp(-i2itvt'4-i2nf »x)  Wab(t',f')|2  - 

"  Xfla ( *jx)  Xbb(*iv,»ix)  .  (106) 

This  is  a  generalization  of  (91). 
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APPLICATION  TO  HERMITE  FUNCTIONS 

This  material  is  heavily  based  on  [5;  appendix  A,  (A-36)  and 

the  sequel].  Let  t  (t)  be  the  n-th  orthonormal  Hermite  function 

n 

with  linear  frequency-modulation,  as  given  in  [5;  (A-36)].  Also 
let  waveforms 


a(  t) 

-  Ck(^t), 

b(t)  ■ 

-  C1(rt),  c(t)  -  Cm(*ut), 

d(t)  -  Cn(rt). 

(107) 

The 

particular 

cross 

WDFs 

wab<t,f) 

-  JdT 

exp( -i2 Jif t)  Ck(//t+Jj/WT) 

C^(Yt-*iYt)  , 

wcd(t,£) 

-  JdT 

exp(  -i2  rtfx)  Cm(/ut+»*/uT ) 

C*(Yt-»5YT)  , 

(108) 

cannot  be  expressed 

in  closed  form.  However 

,  the  cross  WDFs 

wac(t,£) 

■  JdT 

exp( -i2nf x)  tk(//t+»j/ux) 

C*  (  jUt-^/WT  )  - 

«  »  Wkm(/t/t#f//u)  (109) 

and 

Wfed<t,f )  -  ^  Wln(rt,f/Y)  (110) 

can  be  simply  expressed,  in  the  notation  of  [5;  (A-40)  and 
( A— 4 1 ) ] .  Thus,  the  very  complicated  two-dimensional  convolution 
and  Fourier  transform  in  (81),  of  and  WC(J,  can  be  written  in 
a  closed  form  involving  the  product  of  two  generalized  Laguerre 
functions.  Numerous  specializations  are  possible. 
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SUMMARY 

Some  very  general  two-dimensional  Fourier  transforms  of 
convolution  and  correlation  form  have  been  derived  for  various 
combinations  of  WDFs  and  CAFs.  In  particular,  closed  forms  for 
the  convolution  form  are  given  in  (81),  while  results  for  the 
correlation  form  are  given  in  (84).  Numerous  special  cases  may 
be  obtained  from  these  results,  of  which  a  brief  list  has  been 
presented  in  (90)  -  (106). 

Some  extensions  to  more  general  arguments  have  been  derived 
in  appendices  A  and  B.  In  particular,  appendix  A  treats  the  case 
where  a  product  of  CAFs  is  of  interest,  while  the  case  of  a 
product  of  WDFs  is  considered  in  appendix  B.  The  possibility  of 
a  combined  convolution  and  correlation  has  also  been  considered 
in  appendix  A. 

For  signals  reflected  off  moving  targets,  it  is  necessary  to 
define  a  generalized  WDF ,  allowing  for  contracted  arguments. 

This  possibility  has  been  considered  in  appendix  C,  where  a 
two-dimensional  Fourier  transform  and  convolution  has  been 
evaluated  in  terms  of  the  generalized  WDF. 

The  results  of  this  report  should  enable  rapid  evaluation  of 
integrals  of  products  of  WDFs  and/or  CAFs  with  a  wide  variety  of 
arguments  and  including  exponential  terms  with  linear  arguments. 
They  also  significantly  extend  a  number  of  special  cases  already 
known  in  the  literature. 
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APPENDIX  A  -  PRODUCTS  OF  CAFS 

In  this  appendix,  we  will  further  generalize  the  results  in 
(81)  and  (84),  for  products  of  two  CAFs,  to  allow  for  more 
general  arguments.  However,  we  begin  by  considering  general 
two-dimensional  functions  as  in  figure  1.  In  particular,  let 

g(x)  -  xi(va'x)  '  "  x2*vb'r)  '  ( A— 1 ) 

in  ( 4 ) .  Then 

G(f)  »  VV*)  '  H(f)  “  *2(vb'f)  '  (A-2) 

giving 

Jdx'  exp(-i2nfx' )  xx ( vg ,  0r' +ax)  X^(  yr’+ux)  -  exp  (i2nfxttj--"^)  x 
x  Jdf  exp(i2i.f'-c<«r-B„))  *1('-a.r(f'+  j|^))  ^(v6(£'-  2I7))  - 

■  i«tW  exp(12llt'ia2i^)  Jdf'  exP(i2,,£'T)  x 

x  -^4-  +  tI]  *?(vh,  7^1  •  (A-3) 

U  a  cty-0/y  2&J  21  b  ay-0*/  2y) 

Now,  let  va-0v'+a\>,  v'+^v,  where  the  boldface  constants 

are  unrelated  to  their  counterparts;  that  is,  0  need  not  equal  0, 
with  the  same  true  of  a,//,y.  Then  Fourier  transform  (A-3)  on  v' 
to  obtain 
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ll 


dv‘dt'  exp( +i2nv' t-i2nfT' )  (  0v'+otv,|3T'+aT)  x?(  yv'+#iv,  yt'+//t)' 


"p— 1  exp[i2nf'caJi^i]  JJdv'df'  exp(+i2rr\>' t+i2nf 'x)  x 


x  * 


i(^' 


■♦■ov , 


.It 


ar-8/w 


_f 

28 


)  •;( 


yv '  +fi\  r 


ML-  _ 


ay-0/ii 


_f> 

2y, 


( A- 4 ) 


In  general,  we  cannot  proceed  any  further  on  this  double  integral 
of  a  product  of  general  two-dimensional  functions  and  X2* 

Now  let  R^  and  R2  TCFs;  that  is, 

R1(t,x)  -  a(t+J*T)  b*(t-«jx)  -  Rab(t,x)  , 

R2(t,r)  =  c(t+^T)  d*  { t-*jx)  •  Rcd(t,x)  •  (A-5) 

Then  *1  and  *2  become  SCFs: 

♦  l(v,f)  -  tab(v,f)  -  A(f+*jv)  B*(f-4v)  , 

♦  2(v,f)  -  *cd*v'f*  "  c<f+,*v)  D*(f-»jv)  .  ( A-6 ) 

As  a  first  case,  let  y-|3  and  y-0.  Then  (A-4)  becomes 


JJd\>'dx'exp(  i2nv'  t-i2nf  x'  )  Xafa(  0v'+av,  |3t'+cct)  x*d(  0vr+/tiv,  8x'+//x)' 


|  )  |  exP  (j-2nf  T^-^j  JJdv'df'  exp( +i2nv' t+i2nf ' x)  x 

*  +3?|  +‘*ev'+‘sav]  B*(^  +«i|  x 

X  C*(a-"^  "**1  D(^  “*«f  - 
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»  jSP)""1  exp  ^+i2rtf  x5^  -i2nv  « 

x  *ac(f  +l*''<  «-">■!  *bd(f  -*>'*<«--")  -HT(«-„>).  ,A-7) 

Thus,  this  very  general  two-dimensional  correlation  and  Fourier 
transform  of  cross  CAFs  can  be  expressed  as  a  product  of  two 
different  cross  CAFs.  For  8-8-1,  a-a-ij,  this  result 

reduces  to  ( 84 ) . 

As  a  second  case,  let  y»-0  and  y“-0*  Then  ( A—  4  >  becomes 
JJdv'  d-c '  exp(  i2nv'  t-i2nf  x '  )  Xab(av+0v'  ,ax+8x'  )  X*d(//v-8v'  ,/ut-Pt'  )- 

*  JgrwTT  exP  JJdv'df'  exP(  +i2nv'  t+i2nf '  x)  x 

*  A(^  +*s|  -40v'-«4av)  x 

x  c*(^  +H|  -5«^'+^)  D( +,s|  - 

«  |88|-1  exp(+i2nfx^~  -i2nvt-^j  x 

X  ~ac(i  +*it(  «+*» )  ,j|  +JiV{a+/u)j  Wbd  -»}X(  a+,u )  ,^  -^v(a+/w)},  (A-8) 

where  we  used  (61).  Thus,  this  very  general  two-dimensional 
convolution  and  Fourier  transform  of  cross  CAFs  can  be  expressed 
as  a  product  of  two  different  cross  WDFs.  For  a-a-1, 

//-<#- 1,  this  result  reduces  to  (81). 
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As  a  third  case,  let  y-0,  t— 3-  There  follows  a  two- 
dimensional  relation  involving  both  a  convolution  and  a 
correlation: 

jjdv'dx'exp*  i2nv'  t-i2nfr'  )  xab<  @v'+«v,  Pt'+ot)  X*dt -0v'+//v, 

-  |0(3|*  exp  2 nf -i2nvt^jjj  x 

x  WadA[|  +»iv(a+/*)}  W*c^(|  -JjT(a-^/),|  -»jV(«+/#)),  (A-9) 

where  W{t,f)  -  %W(^t,»jf)  again.  Observe  the  conjugates  on 
subscripts  d  and  c  of  the  scaled  WDFs  W. 

For  0-0-1!,  fj-p-l,  this  relation  becomes 

JJdv'dx'  exp(  i2nv' t~i2nf r'  )  Xafe(  v+%v' ,  r+Hr'  )  X*d(  v-Hv' , x+«jx'  >  - 
-  4  exp(  i4nf x)  Wa<J,  (  2t ,  2f+v)  W*c*  (  2t ,  2f-v)  - 

-  exp(i4nfr)  Wad*(  t,f+»*v)  Wbc*(  t,  f-»jv)  .  (A-10) 
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APPENDIX  B  -  PRODUCTS  OF  WDFs 

In  this  appendix,  we  will  also  generalize  the  results  in  (81) 
and  (84),  but  now  for  products  of  two  WDFs,  to  allow  for  more 
general  arguments.  Again,  we  begin  by  considering  general  two- 
dimensional  functions  as  in  figure  1.  In  particular,  let 

g(t)  -  W1(t,f#)  ,  h(t)  -  W2(t,fb)  ,  (B-l) 

in  (4).  Then 

G{ v )  -  *i<v'fa)  '  H(v)  “  *2{v,fb)  »  (B-2) 

giving 

Jdt'  exp(-i2nvt'  )  wx ( pt ' +at , f a )  W*  (  yt '  +vt ,  f  )  -  exp  [i  2  nvta-|-^)  x 

x  Jdv'  exp(i2nv't(  ay-p/v))  *1[yv'+  ^  ,fa]  *2[pv'~  ,fb).  (B-3) 

Now,  let  f a=0f ' +af ,  fb-rf'+^£f  where  the  boldface  constants 
are  unrelated  to  their  counterparts;  that  is,  p  need  not  equal  P, 
with  the  same  true  of  a Then  Fourier  transform  (B-3)  on  f', 
to  obtain 


JIdt’d£ 


'  exp( -i2rtvt'+i2nf ' t)  Wx  (  pt '  +ot ,  pf '  +of  )  W*  ( yt ' +/ut ,  yf 
exp(i2nvt2J±&lj  JJdVdf  '  exp  ^+i2nv' t  (  ay-P/u  )  +  i2 nf  '  xj  x 

x  *1[yv'+  •—  ,pf'+of]  **(pv'-  ,yf'+//f]  .  (B-4) 


41 


TR  8759 


In  general,  we  cannot  proceed  any  further  on  this  double  integral 
of  a  product  of  general  two-dimensional  functions  and  W2 . 

Now  let  R^  and  R2  be  TCFs;  that  is, 

Rx(t,r)  -  a{  t+Jjt)  b*  ( t-Hx )  -  Rab(t,x)  , 

R2(t,x)  -  c  ( t+Jj  x )  d*(t-Hx)  -  Rcd*t,T*  •  (B-5) 

Then  and  *2  become  SCFs: 

♦  ^(Vjf)  -  *ab^v'f*  "  A(f+*jv)  B*(f-%v)  , 

*2(v,f)  -  *ccj^v'f)  *  C  (  f +*{v )  D*  (  f-Jjv )  .  ( B-6  ) 

Substitution  in  (B-4)  yields 

JJdt  'df '  exp(  -  i2nvt '  +  i2nf '  x)  Wafa(  Bt '  +at ,  Pf '  +af  )  W*d (  yt  ’  +/jt ,  yf '  +/if  ) 
-  exp  ^i2nvt^^^j  JJdv'df'  exp  ^+i2xvr  t  (  ay-P//)+i2nf '  xj  x 

x  A^pf '+of+!frv'+l*v/pj  B*  ^0f  ' +af-Jjyv' x 
X  C  ^yf '+|if+4Pv'-J<y/Yj  D  ^yf  ' +/#f-*{Pv ' +  **v/yj  .  (B-7) 

As  a  first  case,  let  y-B  and  y»p.  Then  (B-7)  becomes 
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JJdt'df 'exp(-i2n\»t'+i2nf  'r)  wab(  (St»+«t,  '+«*£)  W*d(  |3t'+/ut,0f '+*/f ) 

-  exp^nvt^^j  JJdv'df'  exp  ^+i2nv'  t$(  a-// )  <-i2nf  •  xj  x 

x  A  ^0f '  +  af +m3v’  +J*v/pj  B*  '  +«f-*s0v' -%v/pj  x 
x  C*  '+//f-Hj(3v'-i<v/|3)  D  '■*■!< v/ej  » 

-  } 03 | _1  exp  (+i2,lvt2lB  "i2BfTS|p)  * 

x  xac(fU_//,+  2I  't(a-")+  ii)  xbd(f(#l“^)“  2I  l|)  *  (B"8) 

Thus,  the  very  general  two-dimensional  correlation  and  Fourier 
transform  of  cross  WDFs  can  be  expressed  as  a  product  of  two 
different  cross  CAFs.  For  8-0-1,  a-a-ij,  *<-//—%,  this  result 
reduces  to  (84). 

As  a  second  case,  let  y--0  and  y--0.  Then  (B-7)  becomes 
JJat-df  'exp(-i2nvt'+i2nf 't)  Wafa{  at+0t ' ,  «f +0f '  )  w*d(//t-pt '  , tit-fit  '  ) 

-  exp^i2nvt^^j  JJdv'df'  exp  (-i2nv'  t|3(  a+/u  )+i2nf '  xj  x 

x  A^0f '+af-Jj0\>'+Vv/8)  B*  ^0f ' +af+Ji$vr -%v/pj  x 
x  C*  (-0f  '-Mf+JiPv'+^v/p)  D  [-0f  '+//f-»i8v'-%v/8}  - 
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■  |  (38 1  ^  exp  (+i2nvt— ^  -i2itf  x 

*  -ac(t*#+^)+  2*  'f(a+'/)4'  18)  W*d(t(a+//)-  jj  .*(*+»*)-  2%)'  <B'9> 

using  (61).  Thus,  the  very  general  two-dimensional  convolution 

and  Fourier  transform  of  cross  WDFs  can  be  expressed  as  a  product  * 

of  two  different  cross  WDFs.  For  p^p-1?,  a-a«l,  //-/#»  1,  this 

m 

result  reduces  to  (81). 

For  t  -  0,  v  -  0,  b(t)  -  a(t),  d(t)  -  c(t),  (B-9)  reduces  to 


JJdt'df '  Waa(at+|3t’  ,af+pf  '  )  Wcc  (  p  t-pt '  ,  //f-pf '  ) 


-  |PP|  1  wac^t(a+//),  f(a+/#)j|2  , 


( B-10 ) 


which  is  nonnegative  for  all  parameter  values  and  waveforms  a(t) 
and  c(t).  This  is  a  generalization  of  (96). 
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APPENDIX  C  -  A  GENERALIZED  WDF 

When  a  signal  is  reflected  from  a  moving  target,  the  effect 
is  to  contract  (or  expand)  the  time  scale  of  the  echo,  rather 
than  cause  a  frequency  shift.  This  requires  us  to  consider  a 
more  general  version  of  a  WDF.  To  begin,  if  waveforms 

a(t)  a  a(ot)  ,  b(t)  a  b(at)  ,  a  >  0  ,  (C-l) 

then  their  cross  WDF  is 

Wab(t'f)  "  a  wab(ott'f/<x)  •  (C-2) 

Thus,  we  have  need  to  consider  integrals  of  the  form 

KaJJdt'df'  exp(-i2nvt'+i2nf 'x)  Wab<t',f')  W*d(t-at',f-f'/a).(C-3) 
This  form  is  general  enough  to  accommodate  integrand 

Wafa( 0t ' , 0f ' )  W*d( t-at ' ,f-f'/a)  ( C—  4 ) 

by  a  change  of  variable. 

To  accomplish  evaluation  of  (C-3),  we  must  define  a 
generalized  WDF  as 

Wab(t,f;p)  a  Jdx  exp(-i2nfx)  a(t+px)  b*(t-(l-p)x)  .  ( C—  5 ) 

Then  we  have  the  usual  WDF  as  a  special  case,  namely 

Wab(t,f;‘«)  -  Wab(t,f)  .  ( C-6 ) 
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Also,  (C-5)  enables  us  to  evaluate  the  following  more  general 
integral  according  to 


J' 


dt*  exp( -i2nf t ' )  a(t')  b  <t-at») 


p  exp( -i2nf tp)  Wab(pt,pf;p)  ;  p  - 


( C-7  ) 


Now  we  are  in  a  position  to  reconsider  integral  K  defined 
above  in  (C-3): 


K  -  JJdt'df'  exp(-i2nvt'+i2nf 'x)  Jdu  exp(-i2nf'u)  a(t'+*ju)  x 
x  b*(tf-»iu)  Jdv  exp{  i2n(  f-f  '/«)v]  c* ( t.~at '  +%v )  d(  t-at ' ~*jv )  .  (C- 


8) 


The  integral  on  f  yields  5(x-u-v/a).  Integration  on  u  then 
yields 

K  -  JJdt'dv  exp( -i2rivt '  +i2nfv )  a( t'+^x-^v/a)  x 
x  b*(  t'-Jjx+Hv/a)  c*(t-at'+*jv)  d(t-at'-^v)  .  (C-9) 


Now  let 


x  -  t'+Jjx-ljv/a  ,  y  -  t'-^x+Jj v/a  ; 

t'  -  *{(x+y)  ,  v  -  a( y-x+x )  .  (C-10) 

The  Jacobian  of  this  two-dimensional  transformation  is  a,  leading 
to 
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X  «  a  JJdx  dy  exp{-iRv(x+y)+i2nfa(y-x+T) ]  x 


if  if 

x  a(x)  b  (y)  c  (t+*jaT-ax)  d(t-haT-ay)  - 


a  exp(i2rcafx)  Jdx  exp[ -i2n(  af+J*\>)x )  a(x)  c*  ( t+lsaT-ax )  x 
x  Jdy  exp[  i2n(  af~^v)y  ]  b*(y)  d( t-^ax-ay )  * 


( 1  +  a) 


expfiZR^^)  W„  1  1  x 

2  1  +  a  )  ac  l  1  +  a  1  +  a  1  +  a; 


X  "bdl  1  +  a  '  1  +  a  '1  +  aJ  ' 


(C-ll) 


by  use  of  (C-7).  For  a  ■  1,  this  reduces  to  alternative  form 
(83),  upon  use  of  (C-6)  and  (61). 


* 
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